"PT-symmetry breaking and universal chirality in a "PT-symmetric ring 
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We investigate the properties of an iV-site tight-binding lattice with periodic boundary condition 
(PBC) in the presence of a pair of gain and loss impurities ±ry, and two tunneling amplitudes 
tojtb that are constant along the two paths that connect them. We show that the parity and 
time-reversal (PT)-symmetric phase of the lattice with PBC is robust, insensitive to the distance 
between the impurities, and that the critical impurity strength for PT-symmetry breaking is given 
by jpt = | to —tb\> We study the time-evolution of a typical wave packet, initially localized on a single 
site, across the "PT-symmetric phase boundary. We find that it acquires chirality with increasing 7, 
and the chirality reaches a universal maximum value at the threshold, 7 = 7pt, irrespective of the 
initial location of the wave packet or the lattice parameters. Our results imply that PT-symmetry 
breaking on a lattice with PBC has consequences that have no counterpart in open chains. 



Introduction. Recently, lattice models with a non- 
Hermitian Hamiltonian that is parity and time-reversal 
(PT)-symmetric, and their experimental realizations in 
optical waveguides [TJ [2] and coupled resistor-inductor- 
capacitor (RLC) circuits with balanced loss and gain [3], 
have become a focal point of research. Properties of 
continuum, non-Hermitian, PT-symmetric Hamiltonians 
were first investigated by Bender and co-workers more 
than a decade ago [4, 5 ; they showed that although 
such a Hamiltonian is not Hermitian, it has purely real 
eigenvalues over a range of parameters. Since then, 
significant theoretical progress on non-Hermitian, PT- 
symmetric Hamiltonians in continuum theories and lat- 
tice models has occurred [SHIP]. Traditionally, the re- 
gion of parameters where the spectrum of the Hamil- 
tonian is real and its eigenfunctions are simultaneous 
eigenfunctions of the combined PT-operation is defined 
as the PT-symmetric region, and the emergence of first 
pair of complex (conjugate) eigenvalues is defined as 
VT symmetry breaking. During the past two years, 
one-dimensional, tight-binding lattice models with PT- 
symmetric, position-dependent tunneling amplitudes and 
on-site potentials have been extensively investigated [TT1 - 
[T4] . These investigations have shown that although 
the PT-symmetric phase is exponentially fragile in a 
constant-tunneling-model |5], it can be strengthened by 
choosing appropriate tunneling profiles [15J and that the 
degree of PT-symmetry breaking in such systems can 
be tuned by the tunneling profile [16 . The equivalence 
between a PT-symmetric, non-Hermitian system and a 
Hermitian system, too, has been demonstrated in limited 
cases [12 El [TBI. 



In the past decade, evanescent ly coupled optical waveg- 
uides [19 have provided a unique realization of a one di- 
mensional lattice with tunable tunneling amplitudes [20 , 
disorder [21 , and non-Hermitian, on-site, impurity po- 
tentials [22]. Today, these systems are a promising can- 
didate for experimental investigations of PT-symmetry 
breaking phenomena in lattice models [TJ [2]. In these 
systems, the tunneling is controlled by the width of 
the barrier between adjacent waveguides and the com- 
plex, on-site potential is determined by the local refrac- 



tive index. The number of waveguides in a typical ar- 
ray is N ~ 10 — 100, and the experimental realizations 
correspond to a lattice with open boundary condition. 
Therefore, most theoretical investigations have focused 
on truncated open lattices, but properties of chains with 
periodic boundary conditions (PBC) have not been ex- 
plored; for a notable exception, see [23]. We remind the 
reader that although the differences between properties 
of open chains and lattices with PBC vanish as N — >• 00, 
for small N < 100, relevant here, they can be significant. 

In this paper, we investigate an TV-site lattice with 
PBC in the presence of a pair of non-Hermitian im- 
purities ±27. Note that on a lattice with PBC, such 
a pair represents PT-symmetric impurities irrespective 
of the individual impurity locations; in other words, for 
given impurity locations, a "parity" operator can be de- 
fined such that the impurity Hamiltonian becomes PT- 
symmetric. The lattice is characterized by two tunneling 
amplitudes that are uniform along the two paths that 
connect the impurities, but may be different from each 
other. 

Our salient results are as follows: i) The PT- 
symmetric phase is robust [24 and insensitive to the dis- 
tance between the non-Hermitian impurities, ii) The crit- 
ical impurity strength jp T at which the PT-symmetry 
breaks is equal to the difference between the tunneling 
amplitudes along the two paths, and is thus widely tun- 
able, iii) As 7 increases from zero, a time-evolved wave 
packet acquires chirality that is independent of the ini- 
tial wave packet location or the distance between the 
impurities, iv) This chirality, quantified by a steady- 
state, dimensionless momentum ^(7) reaches a univer- 
sal maximum at the PT-breaking point, p(jpt) = 1- 
For 7 > 7pt> although the norm of the wave function, 
or equivalently, the net intensity increases exponentially 
with time [25], the chirality is reduced. Our results pre- 
dict that PT-symmetry breaking in a lattice system with 
PBC is accompanied by novel signatures that have no 
counterpart in open chains. 

Tight-binding model We start with a one-dimensional 
chain with N lattice sites and periodic boundary con- 
ditions. Without loss of generality, we take the gain 



and loss impurities (27,-27) at positions (l,d) where 
2 < d < N. The Hamiltonian for this chain is given by 
J^pt = ^0 + V where the Hermitian tunneling Hamilto- 
nian Hq is given by 



N 

H = - ^t(i) (a\ +1 ai + aja i+ ij , 



*(*) 



{ 



i b > 1 < i < d, 

t > d<i<N. 



(1) 
(2) 



Here t(j) is the tunneling amplitude between adjacent 
sites j and j + 1, a] (a^) is the creation (annihilation) 
operator for a single-particle state \j) localized at site j, 
and periodic boundary conditions imply that a N +i = cl[- 
The PT-symmetric, non-Hermitian potential is given by 



V 



ij 



\a\cL! -a\a,d) ^ V ] . 



(3) 



When 7 = 0, the energy spectrum of the Hamiltonian 
Hpt is given by E(k,k r ) = — 2to cos(&) = — 2t^ cos(/c / ), 
and is bounded by 2max(to,t5). When to > t&, the ^ 
eigenmomenta correspond to purely real /c values and £/ 
values that are either real or purely imaginary; when to < 
4, the situation is reversed [26] . When 7 7^ 0, since the 
tunneling amplitudes along the two paths between the 
source and sink are constant, an arbitrary eigenfunction 
|^) = ^2- fj\j) with energy E(k,k') can be expressed 
using the Bethe ansatz as [TTJ [26] 



In 



As'm(k'n) 
Psm(kn) 



- Bcos(k f n) 

- Q cos (kn) 



d 



l<n<d, 
+ l<n< N. 



(4) 



By considering the eigenvalue equation Hpp^) = E\tp) 
near impurity locations, we find that the dimensionless 
quasimomentum k (or equivalently k') obeys the follow- 
ing equation [26 

M (Jfc, k') = tl sin[fc'(d - 1)] sin[/c(7V - d - 1)] 
+ t\ sin[fc'(d + 1)] sin[/c(7V - d + 1)] 
- 2£ 6 £ (sin(fc'd) sin[k(N - d)] + sin(fc') sin(fc)) 
+ 7 2 sin[fc'(d - 1)] sin[k(N - d + 1)] = 0. (5) 

When to > *6? the real quasimomenta k are constrained 
by — 7r < /c < 7r, and states with quasimomenta /c and 
— fc are degenerate except when k = 0,7r, provided the 
corresponding k' is purely real. If fco is a real quasimo- 
mentum, then 7r — fco is also a quasimomentum if and 
only if N is even; when to < *6, the situation is reversed. 
Thus, in contrast with an open chain, the spectrum of a 
chain with PBC has a particle-hole symmetric spectrum 
if and only if N is even. 

Figure [l] shows the numerically obtained typical phase 
diagram for the critical impurity strength jpp as a 
function of fractional distance between impurities /i = 
(d — 1)/N. Note that since there are two paths from the 
source 27 to the sink —27, we restrict the fractional dis- 
tance to < \i < 1/2. These results are obtained for 
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FIG. 1. (color online) Typical VT phase diagram of a lattice 
with PBC. These results are for a chain with N — 40, to = 1 
and impurities ±ry at sites (l,d) for different values of the 
sink-position d and different tb < to; we get identical re- 
sults for odd N or tb > h. Remarkably, the critical impurity 
strength jpt(i^) = \to — tb\ is independent of the fractional 
distance \i between the impurities. Thus, the PT-symmetric 
region in a chain with PBC is more tunable that its open 
chain counterpart. 



a chain with N = 40, to = 1 and t& < to- We obtain 
identical results for odd TV and different values of tun- 
neling amplitudes including tb > to, when the impurity 
strength 7 is measured in the units of max(£o,4). We 
also find that the weak /i-dependence of the critical im- 
purity strength jpp/ max(to, h) vanishes as N increases. 

This remarkable phase diagram predicts that in a chain 
with PBC, the critical impurity strength 7pt(a0 ^ s inde- 
pendent of the inter-impurity distance n, and is given by 
Ipt = \h ~ to\. It implies that the fragile PT-symmetric 
phase in an open chain with constant tunneling [9J [14] is 
stabilized and strengthened by periodic boundary condi- 
tion, and that the critical impurity strength jpp can be 
easily tuned by an appropriate choice of tunneling am- 
plitudes. 

To gain insight into the insensitivity of 7pt(a0 to the 
inter-impurity distance /i for large N ^> 1, let us consider 
Eq.([5} in the limit l«rf«iV, 



sin(/c'/iA0 sin [fc(l - /n)N] 
sin(fc') sin(fc) 



2toh 



[(to-t b ) 2 +7 2 



(6) 



Graphical solutions of Eq.(|6| show that as 7 increases, 
irrespective of /i, two adjacent quasimomenta near k ~ 
7r/2 become degenerate and then complex, leading to the 
PT-symmetry breaking. 

Wave function dynamics. Now we consider the real- 
and reciprocal-space time evolution of a wave function 
that is initially localized on a single site. In an optical- 
waveguide realization of a "PT-symmetric system, this 
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FIG. 2. (color online) The left-hand column shows the evolution of real-space intensity /r(j, £) with increasing 7 > in an 
N — 32 lattice with impurities ±ry at sites 1 and 16, shown by gray circles in the top panel, to = 0.5, £& = 1.0, and initial wave 
packet localized at site mo = 8. When 7 = the wave packet diffuses without preferential chirality. As 7 increases towards 
7pt (center panel) and beyond (bottom panel), the wave function evolution becomes chiral. Recall that due to PBC, the 
bottom-most lattice site in each panel is connected with the top-most lattice site. The right-hand column shows corresponding 
reciprocal-space intensity Im{u^€). The vertical axis, 1 < u < N corresponds to reciprocal space- index — tv < p u < ty. As 7 
increases, we see that the reciprocal-space intensity develops a clear peak at a finite, positive value of p u , consistent with the 
chiral motion of the wave packet in the left-hand column. 



initial state is most easily achievable. For an arbitrary, 
normalized initial state |^(0)), the wave function at time 
t is given by \ip(t)) — ex.-p(—iHpTt/h)\ip(0)} where h is 
the scaled Planck's constant and the time evolution oper- 
ator ex.-p(—iHpTt/h) is not unitary since the Hamiltonian 
Hpt is not Hermitian. We denote the site- and time- 
dependent real-space intensity by Ir(j, t) = \{j\^{t))\ 2 
where j = l,...,iV denotes the site index, and use 
Im(u^) = \(u\ip(t)}\ 2 to denote the reciprocal-space in- 
tensity where the discrete index u = 1,...,7V corre- 
sponds to reciprocal-space index p u = tt(2u/N — 1) with 

-7T <Pu < 7T. 

The left-hand column in Fig. [2] shows the typical evo- 
lution of real-space intensity with increasing impurity 
strength. These results are for a lattice with N = 32, to = 
0.5,4 = 1-0, the source and sink at sites 1 and d = 16 
respectively, and the initial wave packet localized on site 
mo = 8. The vertical axis in each panel indicates the site 
index, and the horizontal axis denotes time measured in 
units of 27rS/max(to,tb). When 7 = (top panel), the 
wave packet diffuses, suffering a change in speed at the 
impurity locations consistent with the change in the tun- 
neling amplitude. As 7 increases towards jpp (center 
panel) and beyond (bottom panel), the wave packet evo- 
lution acquires chirality and the overall intensity also in- 
creases from its 7 = value. In an open chain, there 
is only one path from the source to the sink; in con- 
trast, a chain with PBC has two such paths. Physically, 



the chirality implies that, on average, the path with the 
higher tunneling amplitude is preferred over the other. 
Note that chirality does not represent a preferential flow 
from source to the sink - the wave packet motion contin- 
ues past the sink, to the source again - but rather the 
handedness of the motion. We also emphasize that when 
7 = (top panel), on average, both paths are equally 
preferred. These results are robust, independent of the 
initial location mo of the wave packet or the distance 
(d — 1) between the impurities. 

The right-hand column in Fig. [2] provides a comple- 
mentary view of chirality development, with correspond- 
ing evolution of the reciprocal-space intensity lM(u,t). 
The vertical axis in each panel indicates discrete index 
1 < u < N that translates into p u G (— 7r,7r]. The top 
panel shows that when 7 = 0, as the wave packet dif- 
fuses, its average momentum is zero. As 7 approaches 
7pt (center panel) and beyond (bottom panel), we see 
that the reciprocal-space intensity develops a pronounced 
peak at a finite, positive value, consistent with the chi- 
rality seen in the left-hand column of the figure. In addi- 
tion, results at longer times T ~ ION ^> N show that the 
reciprocal-space intensity distribution for 7 > reaches 
a steady-state. 

To quantify this chirality and to dissociate it from the 
exponentially increasing net intensity Ir(£) = J2 • Ir(j, i) 
that occurs for 7 > jpp [25 , we consider the dimen- 
sionless, discrete momentum operator on the lattice with 



PBC, 



W)\pW)) 



N 

£ 

3=1 



(g*+l+g*)(/j+l-/i) 



(7) 



where \<j>{t)) = E,- <&(*)! J>> |V>(*)> = Ejfj(t)\j), and 
the normalization factor in the denominator is necessary 
due to the non-unitary time evolution. As is hinted by 
the right-hand column in Fig. |2j the momentum expec- 
tation value Pi/j(t) = (i/j{t)\p\ilj{t)) in a given state oscil- 
lates about zero when 7 = and when 7 > 0, reaches a 
steady- state value ^(7) = f Q p{t')dt' jT . Note that since 
IpvWI — 1 f° r an y initial state and time, the magnitude 
of p(j) is bounded by unity. 
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FIG. 3. (color online) Dependence of chirality, quantified by 
the steady-state momentum ^(7), as a function of impurity 
strength 7 for different locations d of the loss impurity —27; 
the gain impurity ij is located on site 1. The initial wave 
function is localized on site mo = 10. The momentum ^(7) 
varies linearly with 7 at small 7/7PT <C 1. It reaches a uni- 
versal, maximum value, p = 1, at the ^T-symmetry breaking 
threshold. For 7/7PT > 1, the steady- state momentum de- 
creases, although the net intensity increases with time. These 
results are independent of the lattice parameters and the ini- 
tial wave function. 

Figure [3] shows the evolution of the dimensionless, 
steady-state momentum ^(7) across the T^T-symmetry 



breaking threshold for different locations of the loss im- 
purity. These results are for a lattice with TV = 32, to = 
0.5, £5 = 1, and thus, 7px/ m ax(£o, h) — 0.5. The initial 
location of the wave packet is mo = 10, and we have used 
normalized time T = 500 to numerically obtain the aver- 
age. These results are generic and qualitatively indepen- 
dent of the initial location tuq or the lattice parameters. 
When 7 = there is no chirality to the wave packet evo- 
lution and p(j) = 0. At small 7, we see that the ^(7) 
increases linearly with 7, as is expected from a first-order 
perturbation theory calculation; the slope of this line de- 
creases monotonically with d. We find that the steady- 
state momentum reaches a universal value, p = 1, at the 
VT -breaking point and decreases linearly for 7 > 7pt- 
We emphasize that although the total intensity increases 
exponentially with time, the chirality, which captures the 
handedness of motion of the wave packet, decreases for 

7 > 7pt- 

Discussion. In this paper, we have investigated the VT 
phase diagram and signatures of T^T-symmetry breaking 
in a lattice with periodic boundary conditions. We have 
presented a model of a lattice with PBC and two uniform 
tunneling amplitudes, and shown that the T^T-symmetric 
region for such a model is robust, insensitive to the dis- 
tance between the loss and gain impurities, and widely 
tunable in size. We have shown that in such a lattice 
with PBC, where there are two different paths from the 
source to the sink, the motion of a wave packet acquires 
a chirality when the impurity strength is nonzero. We 
have predicted that the PT-symmetry breaking in such 
a lattice is signaled by a universal, maximal value for the 
steady-state momentum that quantifies this chirality. 

Traditionally, the investigation of signatures of VT- 
symmetry breaking has focused on dependence of the in- 
tensity profile Iji(j,t) or the net intensity, on time and 
the impurity strength. These quantities vary smoothly 
across the phase boundary [16j [25]. In this paper, we 
have shown that the chirality, on the other hand, shows a 
peak with a universal value at the PT-symmetry break- 
ing threshold. We note that such a measurement will 
require the knowledge of relative phases of the wave func- 
tion values at adjacent sites, ^(7) ex J^ ■ Im(/*/ J+ i), and 
therefore is more complex than the site-dependent inten- 
sity measurements [TJ |2j [I9j Hi] . 

Our results show that T^T-symmetric rings [23] dis- 
play properties that have no counterparts in open chains. 
Their in-depth investigation, including the effects of dis- 
order and non-linearity, will deepen our understanding 
of T^T- symmetric systems with periodic boundary condi- 
tions. 
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